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The sensitivity in optical interferometry is strongly affected by losses during the signal propagation or at 
the detection stage. The optimal quantum states of the probing signals in the presence of loss were recently 
found. However, in many cases of practical interest, their associated accuracy is worse than the one obtainable 
without employing quantum resources (e.g. entanglement and squeezing) but neglecting the detector's loss. 
Here we detail an experiment that can reach the latter even in the presence of imperfect detectors: it employs a 
phase- sensitive amplification of the signals after the phase sensing, before the detection. Since our method uses 
coherent states as input signals, it is a practical technique that can be used for high- sensitivity interferometry 
and, in contrast to the optimal strategies, does not require one to have an exact characterization of the loss 
beforehand. 



From gravitational wave measurements to the realization 
of optical gyroscopes, the estimation of an optical phase (j) 
through interferometric experiments is an ubiquitous tech- 
nique. For each input state of the probe, the maximum accu- 
racy of the process, optimized over all possible measurement 
strategies, is provided by the quantum Fisher information 
H(f) through the Quantum Cramer-Rao (QCR) bound [1, 2]. 
The QCR sets an asymptotically achievable lower bound on 



the mean square error of the estimation 8^ 
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where M is the number of repeated experiments. In the ab- 
sence of noise and when no quantum effects (like entangle- 
ment or squeezing) are exploited in the probe preparation, the 
QCR bound scales as the inverse of the mean photon number, 
the Standard Quantum Limit (SQL). Better performances are 
known to be achievable when using non classical input signals 
13-9]. However, the difficulty involved in implementing them 
implies that all experiments up to now have been performed 
using post-selection, and cannot claim a sub-SQL sensitivity 
when all resources used in the phase sensing are accounted 
for nidi . Additionally, in the presence of loss, the SQL can 
be asymptotically beaten only by a constant factor |[lll - [l3ll . 
so that sophisticated sub- SQL strategies 0, [l4l - [l7ll (imple- 
mented up to now only for few photons) may not be worth the 
effort. This implies that, for practical high- sensitivity inter- 
ferometry, the best resource exploitation (or, equivalently, the 
minimally invasive scenarios) currently entail strategies based 
on the use of a coherent state \a), i.e. a classical signal. Its 
QCR bound takes the form (5^|ql ^ {2Mr]^\a\'^)-^, where 
we consider separately the loss 1 — ^ in the sensing stage and 
the loss 1 — 7^ in the detectors (r] being their quantum effi- 
ciency). Here we present the experimental realization of a 
robust phase estimation protocol that improves the above ac- 
curacy up to ~ (2M^|Qfp)~^, while still using coherent sig- 
nals as input. It achieves the SQL of a system only affected 
by the propagation loss 1 — ^, and not by the detectors' 1 — r]. 

Our scheme employs a simple, conventional interferometric 
phase sensing stage that uses coherent- state probes. These are 



amplified with an optical parametric amplifier (OPA) after the 
interaction with the sample, but before the lossy detectors. No 
post-selection is employed to filter y, llO|] the output signal. 
The 0R\ (an optimal phase-covariant quantum cloning ma- 
chine UM) transfers the properties of the injected state into a 
field with a larger number of particles, robust under losses and 
decoherence [19]. Our approach is suitable for the analysis of 
fragile samples (weak regime of the interaction), since the am- 
plification acts after the interaction of the probe state with the 
sample. A similar scheme was used to perform interferometry 
with single-photon probes ll2Qn . but the small intensity of its 
probing signal could only achieve limited accuracy. 

Theory - The probe is a i^-polarized coherent state 
|<^)i^|0)y, with a = |a|e*^. The interaction with the 
sample induces a phase shift (j) between the polarization 
components and the sample loss 1 — ^ reduces the state 
amplitude to /3 = \^a\ the state evolves to |^^) = 
|e-^^/2^cos((/)/2))i:/|2e-^^/2^sin(0/2))y. The QCR bound 
for this state is ^^^ ^ l/(M2|/3p). In the absence of ampli- 
fication, the detection losses 1 — r] would increase the QCR 
to S'^(l)sQL > l/(M27^|/3p). To prevent this and to attain 
the previous bound, we implemented the operations shown in 
Figdt: a A/4 wave-plate and the 0R\, described by the uni- 



tary UoPA = exp[^(a^^ 
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)/2 + h.c], where g = \g\e 
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is the amplifier gain, and an and ay are the annihilation op- 
erators of the two modes. After these operations, the state 
1^^'^) = ^/opaI^^) is measured by lossy detectors which 
measure the photon number difference D = riH—nv between 
the two modes, with Ux = a^ax. The detector loss 1 — r] can 
be described through a loss map Cr^ acting on 1^^ ), which 
outputs the mixed state p^'^'^. 

The uncertainty on the phase (j) can be evaluated [1] as 
S(j) = cr((^))|-^^|~^, where (D) is the expectation value 

of I^ on p^'^'^. A calculation of the sensitivity S = 6(j)~^ 
of the whole procedure (see Supplementary Material) shows 
that S depends on the value of the phase (j) to be estimated. 
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FIG. 1. Experimental scheme and theoretical yields, (a) Scheme of 
the amplifier based protocol. The two red lines represent the polar- 
ization modes H and V\ U(j) introduces the relative phase on the 
state; C^ and Cr^ are the loss transformations for the sample (which 
transforms the probe state as |a) -^ |/3) = |v^<^)) and detector 
loss respectively; the blue box is a A/4 plate and the orange box 
represents the optical parametric amplifier (OPA); the black and the 
gray boxes represent the detection and data-processing. The quantum 
Fisher information ifampi is achievable optimizing over the detec- 
tion and the data-processing; the Fisher information /ampi is achiev- 
able optimizing over the data-processing; the sensitivity S is what is 
achieved for a given detection and data-processing, (b) Plot of the 
ratio between the sensitivity square (5'max)^ and the quantum Fisher 
information i/sQL connected to the standard quantum limit (SQL) 
as a function of the nonlinear gain g of the amplifier and of the de- 
tection efficiency 77, with |/3|^ = 20. Our scheme achieves the SQL 
for a wide range of parameters, (c) Plot of the ratio between the 
quantum Fisher information ifampi from our scheme and if sql as a 
function of the nonlinear gain g of the amplifier and of the detection 
efficiency 77, with |/3|^ = 20. (d) Comparison between the classical 
Fisher information /ampi (points) and the sensitivity (5'max)^ (lines) 
for |/3|2 = 9. 



The maximum sensitivity is obtained forcj) = 7r/2 by setting 
A = 20, i.e. 



^ V^(l + 2n + 2v/n(l+n)) 



2n(l + 7? + 27^n) + \f3\^ [l + 2n + r]n{6 + 8n)] 

(1) 

where n = sinh^ g. Because of the dependence of S on 0, 
to achieve the maximum sensitivity 6'max an adaptive strategy 
112 ill is necessary. In the Supplementary Material we show 



that it is sufficient to use a simple two- stage strategy in which 
we first find a rough estimate of the phase 0est employing con- 
ventional phase estimation methods, and then we use it to tune 
the zero-reference so that our scheme operates at its optimal 
working point detailed above. We also show that the resources 
employed in the first stage of this adaptive strategy are asymp- 
totically negligible with respect to the resources employed in 
the second high-resolution stage. We will thus neglect the first 
stage in the following analysis. 

Efficiency of the phase estimation - To gauge the efficiency 
of our method we start by noting that for n :^ (2//)"^ and 
|a| > 1/(^2^) the sensitivity ([B gives Smax — V2\f3\, so 
that in this regime the QCR bound (5^0 ^ l/(M2|/3p) of 
the state |^^) (before the amplification and the detector loss) 
can be attained by our detection strategy. More generally, our 
detection strategy is also optimal to extract the information on 
the final state p^ '^ in other regimes. In fact, the quantum 
Fisher information if ampi of this state is (see Supplementary 
Material) 



i^ampl(|/3|,^,7?)=2|/3p7? 
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y^l ^ 4r]{l - r])r 
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where ^eff = l/41og[(7^e^^ + 1 — r])/{r]e~'^^ + 1 - r])], and 
where again we maximized the ^-dependent quantum Fisher 
information by choosing (j) + A/2 — 6 = 7r/2 (implying that 
again an adaptive strategy must be employed). The quantum 
Fisher information, through the QCR bound, measures the 
best precision achievable when optimizing over the possible 
detection and data processing strategies (see Fig[T^). To show 
that our choice to measure the photon-number difference D 
can be optimal, note that for n :^ (87^) ~^ and |/3p :^ 1/2 the 
ratio S'^/Hampi -^ 1, see Fig[Th-c. In addition, our data pro- 
cessing can be optimal for even a wider range of parameters. 
In fact, the sensitivity S closely tracks the classical Fisher in- 
formation /ampi also for small values of n, see FigdJi. The 
quantity /ampi represents the maximum amount of informa- 
tion that can be extracted from the probe state using our choice 
of measurement, optimizing over all possible data-processing. 
We now compare our method to other strategies, using as 
a benchmark the SQL (^^^sql ^ l/(M2|/3p), which would 
be achieved by a probe coherent state with |/3p average pho- 
tons using lossless detectors. As detailed above, our method 
asymptotically (for sufficiently large n and |a|) matches this 
strategy achieving shot-noise limited sensitivity, see dash- 
dotted line in FiglJ^. In other words, increasing the ampli- 
fier gain, the effects of the detector loss can be asymptotically 
removed ll22ll . Consider now the case with no amplification, 
where a coherent state is subject to both the sample and detec- 
tor loss (FiglJt). This is the strategy conventionally used in 
interferometry. Our method clearly always outperforms it, see 
Figsllt-d and the continuous line in FiglJ^. Recently, the op- 
timal strategy in the presence of loss was derived [5] (Fig|2t)). 
It employs the state that maximizes the quantum Fisher infor- 
mation in lossy conditions. Of course, this strategy cannot be 
beaten if one could access the optimal measurement that at- 



tains the QCR bound. Even though elegant proof-of-principle 
experiments exist |9], both this measurement and the cre- 
ation of these states without using post-selection are beyond 
the reach of practical implementations for the foreseeable fu- 
ture, especially for states with large average photon-numbers. 
In addition, the form of these states strongly depends on the 
value of the loss 1 — ^: this may not be available in prac- 
tice and its experimental evaluation typically requires irradi- 
ating the sample which removes the advantage of using the 
optimal minimally-invasive states. In contrast, our amplifier- 
based protocol uses readily available input states and detection 
strategies, and does not require a priori knowledge since the 
choice of the coherent state is independent of the value of the 
loss. Since our method is devised especially to counter the 
detector loss 1 — 7^, we compare the performance of our states 
with the optimal state calculated for the total amount of loss 
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FIG. 2. Comparison with other schemes, (a) Conventional (unam- 
plified) coherent- state interferometry with sample and detection loss 
1 — ^ and \ — r\ respectively: it can achieve the SQL bound connected 
to the quantum Fisher information (QFI) H^^^. (b) Interferometry 
based on the states that optimize the QFI in the presence of loss, pro- 
posed in ill], with the corresponding QFI, H^j^. (c-d) Comparison 
between our scheme and the conventional scheme of FiglJ^: loga- 
rithmic plots of the ratio (Sampi)^ /Hsql as a function of 77 and g 
for |/3|^ = 20 (c), and as a function of |/3| for g = 3.4 (d). (e) Com- 
parison between the QFI for the above strategies for |/3|^ = 20 and 
g = 3.5. Black dotted line: reference, chosen as the shot-noise 
HsQL. Blue dash-dotted line: ratio between the shot-noise ifsQL 
and the QFI of our method, i/ampi- Red solid line: comparison be- 
tween the shot-noise limit ifsQL and the coherent state phase esti- 
mation with loss of FiglJ^, ii^sQL- Green dashed line: ratio between 
the shot-noise Hsql and the optimal strategy of FiglSt, H^ [5^. 



1 — ^rj, showing that our method can achieve better perfor- 
mance for the practically-relevant case of low values of r] (see 
dashed line in Fig|2^), where the detection strategy is clearly 
not optimized to achieve the QCR bound of the optimal states. 
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FIG. 3. Experimental setup for the practical implementation of the 
protocol. The excitation source is a Ti:Sa laser system, consisting 
in a Ti:Sa mode-locked Mira900, whose output beam is injected into 
the Ti:Sa RegA9000 amplifier. The overall laser system can output a 
1.5W beam at wavelength A = 795 nm. In a first nonlinear crystal, 
the output field is doubled in frequency through a second harmonic 
generation (SHG) process to generate the experiment pump beam at 
wavelength Xp — 397.5 nm of power P — 650 W. The remainder 
of the 795 nm beam is then separated from the pump beam through 
a dichroic mirror, and is prepared in the coherent state |a)+ by con- 
trolled attenuation, spectral filtering (IF) and polarizing optics. The 
coherent state probe then acquires the phase shift by interacting with 
the sample (in our case, a Babinet- Soleil compensator), and is then 
injected into the OPA after the acquisition of the phase. 

Experiment - We now describe the experimental implemen- 
tation in highly lossy conditions, showing that we can achieve 
a significative phase- sensitivity enhancement with respect to 
the coherent probe based strategy. The optical setup is re- 
ported in Fig. [3l To acquire the phase shift to be measured, 
the probe coherent state is injected into the sample, which is 
simulated by a Babinet- Soleil compensator that introduces a 
tunable phase shift between the H and V polarizations. Sub- 
sequently, the probe state is superimposed spatially and tem- 
porally with a pump and injected into the OPA. In this ex- 
perimental realization the phases of the pump and of the co- 
herent state are not stabilized: this will reduce the achievable 
enhancement by a fixed numerical factor of 4. In contrast to 
previous realizations of parametric amplification of coherent 
states 1I23I425I1 which focused on the single-photon excitation 
regime, we could achieve a large value for the nonlinear gain, 
up to ^ = 3.3, corresponding to a number of generated pho- 
tons per mode n ^ 180 in spontaneous emission. In addi- 
tion, our scheme is also able to exploit the polarization degree 



of freedom. After the amplification, the two output orthogo- 
nal polarizations were spatially divided and detected by two 
avalanche photodiodes. Their count rates are then subtracted 
to obtain the value of {D), and recorded as a function of the 
phase 0, introduced by the Babinet. 
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of noisy detectors. This approach involves coherent states as 
input signals, thus not requiring any a priori characterization 
of the amount of losses, and phase sensitive amplification af- 
ter the interaction with the sample and before detection losses. 
By adopting the amplification stage, the accuracy of the pro- 
tocol can reach the performances of a lossless probe state even 
in presence of imperfect detectors. Furthermore, we presented 
the experimental implementation of this protocol in highly 
lossy scenario, showing the advantage of this technique even 
in a highly detrimental regime. As a further perspective, the 
present strategy could be exploited in phase estimation proto- 
cols with noisy detectors involving different classes of probe 
states, including quantum resources such as squeezed light. 
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FIG. 4. Experimental results, (a) Fringe patterns for the amplified 
coherent state {g — 3.3) and (b) for the unamplified coherent state 
strategy for |/3|^ ~ 22.8. (c) Fringe pattern visibility and (d) exper- 
imental enhancement ^exp as a function of the nonlinear gain g for 



|/3|2 - 5.8, ry - 1.46 x 10" 
theory: black solid line) and 



(experiment: black diamond points; 

3\^ - 22.8, T] - 3.48 x 10"^ (ex- 



periment: green star points; theory: green dashed line). Inset: ex- 
perimental plot of the sensitivity as a function of the efficiency ratio 
?7/?7o, where 770 is a reference efficiency. [See Supplementary Mate- 
rial for details on the theoretical plots in (c)-(d).] 

The results of the experiment are reported in Fig. IH Since 
the sample losses 1 — ^ act as a scaling factor on the coher- 
ent state amplitude as a ^ /3 = ayT' we evaluated |a| by 
estimating the average number of photons |/3p after the in- 
teraction with the sample. In FigsE^-b we report the fringe 
patterns obtained by measuring the single counts at the single 
photon detectors for the amplified case and the coherent state 
case, respectively. An enhancement of ~ 200 in the counts 
rate for the former case is observed without significantly af- 
fecting the visibility of the fringe pattern (Fig 14]:), leading to 
an increased phase resolution. We measured the enhancement 
£^exp achievable with our protocol with respect to the conven- 
tional unamplified interferometry, defined as the squared ratio 
between the measured sensitivities with and without the am- 
plifier, see Figllli. Our measurement shows a good agreement 
with the theoretical predictions. A significant enhancement 
up to a value of £^exp = 186.3 ± 9.3 has been achieved. In 
the operating regime of our experimental implementation the 
sensitivity of ^exp scales as t^^/^ (see Fig. [l]:, inset), as for 
a coherent state only. Hence, the observed enhancement is 
mainly due to the strong increase in the counts rate due to the 
amplification process. 

Conclusions and perspectives - In this paper we presented a 
strategy to perform phase estimation protocols in the presence 
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In this supplementary material we elaborate on the mate- 
rial presented in the main text, giving more details on the ex- 
perimental procedure and carefully deriving of the formulas 
presented there. In Sec. H we describe the experiment and 
the evolution of the quantum state of the probe as it evolves 
through the apparatus. In SecJIlwe derive the phase sensi- 
tivity S of our apparatus, Eq. (1) of the main text. In Sec. Hill 
we give the details of our simple two-stage adaptive scheme, 
showing how the first stage (where a rough estimate of the 
phase (j) is recovered) can be neglected asymptotically, as it 
requires asymptotically vanishing resources. In Sec. [Ivl we 
calculate the explicit form of the output state of our appara- 
tus. In Sec.|Vlwe calculate the quantum Fisher information of 
the output state, and in Sec.|Vllthe classical Fisher informa- 
tion that results from fixing the detection scheme to the one 
we employ in the experiment. Finally, in Sec. I VIII we give the 
details of the theoretical model we employed to analyze the 
experimental data. 



I. EXPERIMENT 

The scheme of the proposed apparatus is sketched in Figdl 
The probe is a horizontally (H) polarized electromagnetic 
field prepared in the coherent state \a)H\0)v with a = | a | e*^ . 
It is sent through an interferometric setup to interact with the 
sample. The sample induces a phase shift ^ on the system and 
is characterized by a loss 1 — ^. The aim of our apparatus is 
to determine (/), while employing a low intensity signal. The 
phase shift is induced through a unitary transformation of the 
type 



U, 



= e-^(«- 



O^a_0 



(1) 



where a_ = {an — av)/V^ is the annihilation operator con- 
nected to the — polarization. The loss is induced through a 
completely positive map C^ of the form 

l^e-l _ 1 )n/2 (jt(j 

AW = ^AnpiAnY, with An = ^ ^=^a^^—(2) 



n\ 



where p is an arbitrary state. Since the action of the phase 
unitary /70 and of the loss C^ commute, we can consider these 
two as independent processes that occur during the interaction 



with the sample. The action of the loss map on a coherent state 
simply shifts its amplitude >C^[|a)(a|] = \^a){^a\, without 
changing the form of the state. Thus, our choice of coherent 
state probes will not depend on the noise characteristics of the 
sample. Consider first the unitary part of the interaction t/^: 
the state evolves as 



m) = U^\a)H\0)v = 



= |e-'*/2acos((/)/2))i/|ze-'^/2^sin((/)/2))v. 



(3) 



Then, the action of the loss £^ reduces the ampUtude of the 
coherent states so that, after the interaction of the sample, the 
probe has evolved to 



\<) 



'^/2/jcos((/./2))^|ze-'^/2^sin((/./2))v, (4) 



with P = \^a. Before the amplification, a relative phase- 
shift of 7r/2 is inserted between the H and the V polarization 
components by means of a A/4 birefringent waveplate with 
optical axis oriented at 0°, leading to: 

|e-^^/2/3cos((/)/2))^| - e-^^/2/3sin((/)/2))y. (5) 

The resulting state is then injected an the optical parametric 
amplifier (OPA). The interaction Hamiltonian of the OPA is 

y^OPA = ^hx (a^^al_ ] +H.C. = ihx (a^H ~ ^v ) /2 + H.c. 

(6) 
where a± = {aH ± ay)/>/2, and x is the parameter that 
quantifies the strength of the interaction. It corresponds to a 
unitary operation 

^OPA = exp[^(a^^ - a(.^)/2 + h.c] (7) 

where g = \g\e'^^ = x^ is the amplifier gain (t being the 
interaction time). Form the form of the unitary in (|7]), it is 
clear that the OPA is equivalent to two single-mode squeezers 
acting independently on the modes H and V with opposite 
phases, namely /7opa = SH{—g) ^ Sv{g), where Si{g) = 
exp[-gal ^ /2 + h.c], 1 = H,V. 

After the amplification, the state has evolved to |^^'^) = 

[/qpaI^^). Finally, it is detected by lossy detectors, 
parametrized by a quantum efficiency r]. These are equiva- 
lent to perfect detectors that measure the number of photons, 
preceded by a loss map >C^ 1 1]. The action of this map on the 
state 1^^'^) produces the mixed state 



^f3,g,rj _ 



f3,9\/.Tff3,9\ 



^.[i*r)(*ri] • 



(8) 
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The explicit form of this state will be calculated in Sec.HVl 
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FIG. 1. Scheme of the different steps of the ampUfier-based phase 
estimation protocol. 



II. THEORY OF THE SENSITIVITY OF THE PROTOCOL 

In this section we report the details of the calculation of 
the phase sensitivity S for the proposed apparatus. It is con- 
venient to work in the Heisenberg picture. To this end, we 
need to consider the time evolution of the field operators due 
to the OPA and of the loss map Cn . The latter is equivalent to 
the insertion of a beam-splitter of transmissivity r] along the 
transmission path of the field, seeded by the vacuum state in 
the other input port. By combining the resulting equations for 
the time evolution of the amplifier and of the beam-splitter we 
obtain the following expressions for the field operators at the 
detection stage 






V^ (a^^C + e-'^anS) - i^l^h^H (9) 
-'^avS] - iJl^hl (10) 



Vv(a^vC 



where hn and by are the annihilation operators for the second 
input port of the beam-splitter, C = cosh \g\ and 5* = sinh \g\. 
The chosen strategy to extract information on the phase shift 
(j) is to measure the output photon-number difference D = 
cIjCh — CyCv and to extrapolate the value of from the de- 
pendence of {D) on it. By exploiting the expressions ([9lfT0l) 
for the field operators, the average of D on the state p^'^'^ is 



(D) = 7^|ap^[cos0(l + 2n) + cos(0 + A-2l9)2Vn(rT^] 

(11) 
To evaluate the resolution Scj) on the estimated phase according 
to standard estimation theory, we need to calculate the fluctu- 
ations cr{{D)) on the detected signal. Such quantity can be 
evaluated according to a'^{{D)) = {D'^) — {D)'^ . By evaluat- 
ing the average values {{c^jjChY) and ((cyCy)^), we obtain: 

cr^{{D)) = 7? [a(n, 7/) + cos 0cos(0 + A - 20)b{n,r])] (12) 

where: 

a{n, T]) = 2n(l + r? + 2r]n) + |a|^^ [l + 2n + r]n{6 + 8n)] 

(13) 
(14) 



6(n, T]) = 2^/n{l^n)\a\'^^{l + 7? + 4r]n) 



We note that both the signal and the fluctuations depend on the 
phase difference between the coherent beam and the pump 
beam A. Finally, the resolution of this detection strategy can 



be evaluated according to standard estimation theory as 

d{D) 



S4>=Va^m)/ 



(15) 



\/a{n^ rf) + cos (j) cos((/) + A — 20)b{n^ rf) 



lofpy^^l cos(/)(l + 2n) + cos( 



• + A-2(9)2v'n(l + n)| 
(16) 



The sensitivity is defined as ^ampi = ^0~^. Its optimal oper- 
ating point is achieved for A — 2^ = and for a value of the 
actual phase of (j) = 7r/2, corresponding to the steepest point 
of the signal (D). The sensitivity of the scheme in this point 
is then 



^max — 



|q^P^V^(1 + 2n + 2y^n{l + n)) 



(17) 



The fact that the sensitivity 6'ampi depends on the parame- 
ter (j) we want to estimate implies that the maximum sensitivity 
*5'max can be achieved only by employing an adaptive strategy, 
where some initial measurements are performed to get an esti- 
mate of (j) so that the apparatus can be employed in its optimal 
working point around (/) = 7r/2. 



III. ADAPTIVE PROTOCOL 

In this section we detail a simple two- stage adaptive 
scheme, where first a rough estimate of the parameter (j) is 
found, and then this estimate is employed in a second high- 
resolution stage of the protocol. 

Let (j) be the parameter we want to estimate (the phase) and 
assume that it is encoded in two different families of states, 
i.e. the family {p^}^ and the family {cr^}^. For example, 
the first family can be identified with the states of the system 
at the output of the interferometer when no amplification is 
used. The second family instead is identified as the the state 
at the output of the interferometer when the amplifier is ac- 
tive and where we have set the phase reference in such a way 
that the apparatus gives optimal performances for ^ = 0. In 
what follows we will consider a two stage estimation strategy 
in which i) first we perform Mi measurements on the state p^ 
of the first family to get a preliminary estimation of (j), and 
then ii) we perform M2 measurement on the state a(f) of the 
second family to improve our estimation (of course in the sec- 
ond stage we are facilitated by the fact that we have already 
acquired some info on (j)). 

Let then x = {xi^X2j- ■ ■) the data extracted from the first 



set of measurement and 



a(^i) 



(x) the estimation function we 



use to get the preliminary estimation of (j). Using the quantum 
Cramer-Rao (QCR) bound we have 



^2(/>i = ^Pi(x)[(/>- 



\x)]' 



> 



1 



MiHiicjy) ' 



(18) 



where Pi {x) are the probability of getting the outcomes x 



(g)M: 



when measuring p^ ^ and Hi {(j)) is the quantum Fisher info 



associated with the family {p^j^. For the sake of simpHcity 
we assume that x]^^^ (x) is unbiased, i.e. 



^p,{x)[^-4^;\x)]=o, 



(19) 



(generaHzation to the general case are possible). 

In the second stage of the estimation we use the family 
{cr^jc/, where we modify the way the phase is mapped by 

rescaling it by (j)l^^^{x). This is possible for instance by 
changing the initial phase reference which effectively shifts 

the unknown phase ^ to x = ^ — ^Lt (^)- ^^^^ ^^ ^^^ ^^^ 
parameter we wish to recover. In the second stage, we perform 
measurements on a^^^ obtaining the data y = (^i, ^2, • ' ' )• 

We determine x via the estimator xL/ (^) which again we 
assume to be unbiased, i.e. 



yZP2{y)[x-xi^t\y)]=^ 



(20) 



(here P2 (y) is the probability of getting the outcomes y when 
measuring a^^'^). The whole process can be described hence 
by introducing a joint estimator function 

^i7''''\x,y) = <Pi^,^\x)+x'^^\y) . (21) 

characterized by a probability distribution Pi{x)P2{y) and 
which (by construction) is unbiased, i.e. 

J2 Pi {S)P2 {y) ~& ''^'\x, y) = ct>- (22) 

Let us now compute the variance of the error associated with 
such estimator. Formally this is given by 

x,y 

= J2p^{x) J2P2(y)i^-^i'^i''^'\^^y)f 

X \_ y 

= J2P^ix) J2P^(y) [<^-<P'J:ft'\x)-x'J^^\y)f 

X \_ y 



X 

^ypiix) 



1 



M2H2{x) 



Y,p,{x) 



1 



M2H2{ 



&\^)) 



(23) 



where we used the QCR bound on the estimation of x and 
where -ff2(x) is the quantum Fisher info of the state a{x)- 
The above expression can now approximated by using the fact 

that for sufficiently large Mi, 4'ext \^) — 0. i-e. X — 0. This 
allows us to expand H2 (x) around x = 0, i.e. 



H^ 



<P['^l\x)) C. ^2(0) + (0 - (Ai^^H^)) ^2(0) 



(Mi).-jxx2 TTin 



+i<i>-€:rix)rmio)/2 



(24) 



which yields 



1 



M2 



^'^ ^ w E^i(^) V ^2(0) + (0 - 0i^^\f)) Hi,{o) 



-{4>-4>^-:^i\x)?H'mn 



Y,P^{x) 



M2H2{0) 



e^(.))f|| 



-{^-^'j^;\x)r ""^^'^ 



-{^-&m' 



1 



2H2{0) 



n2. 



M2H2{<d) 



H2{0) 
V2J?2(0) 



Hm 



^2(0) J 



where we used Eq. ( fT9l ) and the definition of S'^cpi. Suppose 
now that H2 (x) achieves its maximum for x = (this is what 
happens thanks to our new choice of reference). This implies 
that Hm = and H^'{0) < 0. Therefore we get 



5^;^ 



1 



M2H2{0) 



1 



5^<P: '"' 



> 



1 



g^(0)l 
' 2JT2(0) 

H^m 



M2H2{0) I 2MiHi{4>)H2{0)l ' 



(25) 



where in the last inequality we used the QCR bound dTSl) . 
Defining M = Mi + M2 the total number of measurements, 
we can write 



J20> 



1 



{l-p)MH2{0) 



\H^m\ 



2pMHi{(P)H2{0) 



,(26) 



with p = Ml /M begin the fraction of measurement we em- 
ploy in the first step of the protocol. This equation provides 
the corrections to the accuracy we get when we adopt the 
adaptive strategy. 

Observation I: It is worth comparing the above bound with 
the accuracy one could get if instead of performing the prelim- 
inary step one could have used all M copies to perform only 
the estimation on the states a^. In this case the resulting accu- 
racy would be 1/{MH2{(I))). Do we gain something by going 
true the adaptive result? A positive answer would require 



1 



l^^'(0)| 



{l-p)MH2{0) L 
which can be cast as 



2pMHi{(l))H2{0)i MH2 



^ 



kV) 



P^A 
p{l-p) 



^B 



(28) 



1^2^(0)1 



r . Since by 



with B = H2{0)/H2{4>) and A = 2MH,WH,io)- 
assumption B ^ 1 and A ^ 0, one can easily verify that 
there are value of p which allows one to obtain Eq. (l27l) if B 
is sufficiently large. 



Observation II: For fixed M we can optimize the right- 
hand-side of Eq. ([26l) with respect to p. This yields 



Popt 



V^42+A-A, 



(29) 



(notice that this is and increasing function of A which is 
always positive and smaller than 1/2 - the latter being the 
asymptotic value reached for A >> 1). Consequently we 
can write 



5^^^ 



1 



1 



l^^^(0)| 



> 



{l-p)MH2{0) 

1 
" {l-p)MH2{0) 
1 



2pMHi{^)H2{0) 

A 

+ - 
p 



MH2{0) {^/WTA-A){l^A-^/A^TA) ' 

0. Therefore we can 



Now, for M >> 1 we have that A 
write 



^'(/>^ 



1 



MH2{0) 

1 
MH2{0) 



[1 + 2VA] 




2|i7-(0)| 



MHi{^)H2{0) 



(30) 



This implies that the resources Mi employed in the first stage 
of the protocol can be neglected, and the precision asymptot- 
ically approaches the QCR of the second stage: the term with 
the square root in ([3Ql) is asymptotically negligible. 



IV. STATE EVOLUTION 



The action of the squeezing operators and of the displacement 
operators can be now inverted according to 

D{a)S{g) = S{g)D{a+) (34) 

S{g)D{a) = D{a-)S{g) (35) 

where a± = acoshg ± a*e*^ sinh^. Using Eqs. (l34H35]) we 
can write 



Si{gH)Di{l3i)\0) = Di{ji)Si{gi)\0) 



(36) 



with ji = f3i cosh gi — f3fe^^^ sinh gi. The output state can be 
then written as 



^/3,g,v 



Cr,{DH{lH)Dv{lv)SH{gH)Sv{gvm{0\ 

Sl,{9H)Sl{gv)D^H{jH)Dl{jv)} 



(37) 



By interchanging the action of the loss >C^ and of the displace- 
ment operators Di{'^i), we obtain 



J^9 



^ = DH{lH)Dv{lv)Cr,[SH{gH)Sv{gvm{^\ 

sU9H)SU9H)}D^H{lH)Di{^v) 



(38) 



where ^i = -^/rjji . Finally, by exploiting the identity (IB2b of 
Appendix ini involving the action of £^ on squeezed vacuum 
states, we can express the output state after detection losses in 
the Gaussian form 

P^'"'" = DH{lH)Dv{jv)SH{gf)Sv{gf) [p'^{N,f,)^ 

(39) 
The expressions for gf^ and N^^ are reported in Eqs. (IB3IIB6b . 



In this section we calculate the explicit form of the output 
state p^'^'^ of our scheme, by exploiting some operatorial re- 
lations for Gaussian states. This will be useful to evaluate the 
quantum and classical Fisher informations in the following 
sections. The state impinging at the measurement stage after 
detection losses can be written in the form: 

4'^'^ = Cr^{SH{gH)Sv{gv)C^ [DH{aH)Dv{avm{0\ 

DUo^H)Dl{av)\sl,{gH)Sl{gv)} 

(31) 
where Di{ai) = exp(a^a^^ — a^ai) is the displacement oper- 
ator such that D{a)\0) = \a). The action of the lossy channel 
^ and of the displacement operators can be interchanged as 



C^[DH{aH)Dv{av)\0){0\Dl{aH)Dl{av) 
= DH{f3H)Dv{M\0){0\DUMDUM 
where f3i = \^ai. The output state then reads: 

f//'" = C^{SH{gH)Sv{9v)DH{f3H)Dv{f3vm{0\ 
Dlj{MDUMSU9H)sUgv)} 



(32) 



A. Eigenvalues and Eigenvectors 



From Eq. (|39l ) one can calculate the spectrum of eigenval- 
ues and eigenvectors of p^'^'^. As a first step, we observe that 
the density matrix of the state takes the form of a separable 
state p\ ^ p\ \ where 



.(0 



Di{ji)Si{gf)pf{Nf)SJ{gf)DJ{ji), (40) 



with / = H^V. Since the state for the two modes has the 
same Gaussian form, the joint spectrum can be obtained by 
analyzing directly the p\^ single-mode state. By expanding 
the density matrix in the Fock basis we obtain: 



n=0 ^ ^ ^ 

\nh{n\SJ{gf)DJ{j,) 

The eigenvalues and the eigenvectors of the state 

Z^n ^n \i^n )i{i^n \ are then respectively 



(41) 



.(0 



,(0 - 



(Nfr 



(33) 



(1 + 7Veff)n+l 



(42) 
(43) 



Finally, the eigenvalues and the eigenvectors of the joint two- 
modes density matrix can be written as 



^I3,9,r] 



Y^ gmJ"^ 



m.n ^ m.n/ 



)w(*m,n| (44) 



„ _ JH) JV) 

fm.n — fm en 



(45) 
(46) 



V. QUANTUM FISHER INFORMATION 

In this section we describe the calculation of the quan- 
tum Fisher information (QFI) of the output state p^'^'^ of our 
scheme. 

The QFI for a generic mixed state a = ^^ crm\Cm){Cm\, 
as reviewed in the Appendix lAl in Eq. (IA4b . can be evaluated 
asH: 



Hch 



(^0^p)^ 



n,m 



-n,m\ \Sm|^0Sn/ 1 



(47) 



Here am and |Cm) are respectively the eigenvalues and the 
eigenvectors of the density matrix, and Cn^m = {c^n — 
crm)'^/{crn + cTm)- In the case of the output density matrix 
p^'^'^ of the amplifier-based protocol the eigenvalues and the 
eigenvectors are parametrized by the indices (n^m), and the 
QFI is 



H{a,^,{gi},{Xi},r]) 



p,q=0 



bQp,q) 



Qp.^ 



(48) 






-2 2^ ei,^-^,n|(^i,i|90^m,n)|' 
i,jf,m,n=0 



where 



\Qi,j Qm,n) 



Qhj 



We observe that, for the density matrix 



^I3,9,r] 



the eigenval- 
ues gm,n (l44l) are independent on the phase (j), and hence the 
first term in Eq.(|48l) vanishes. In order to calculate the sec- 
ond term, it is necessary to evaluate the following quantity: 

I {"^ij |90^rn,n) p. Such term can be written as 

(50) 
Since the eigenvectors for the two-modes present an anal- 
ogous form, it is necessary to evaluate only the term 

i{iIj\ \d(f)ilJrn)i- Let us focus on the I^^V^m )^ state vector. 
Since the dependence on (j) of the state is included only in 
the displacement operator Di{ji), we can write: 



\d^^i^)i = [d^Di{^i)]Si{gf)\m) 



(51) 



The latter can be evaluated by differentiating the displacement 
operator written in normally-ordered form: 

d^Diiji)] = d4e--^\^^\\^^^'e-^>^] (52) 

By differentiating the three exponential with respect to 0, and 
by exploiting the following commutation relation: 



the derivative of Di{^i) reads: 



(53) 



(54) 
The scalar cj^^ ^^ ;,, ,^^^ and the operator F^^^^^ ^^^ _^_^ (a; , a] ) 
are respectively: 



C. 



^aig,,A,,,,,,^(«;' 4) = id<t>7i)a] - {d^lt)ai (56) 

By replacing the latter expressions in Eq.dSB. the scalar prod- 
uct ; {'ip\ ' I dcji'ipm ) I can be evaluated as: 



i{i^1\d^i^^:;})i= MsKgf)D\{ii)[cf>^, 



^,^z,Az,77,0 






(57) 

Such average value can be evaluated by exploiting the opera- 
torial identities 

S\g)aS{g) = acosh^ — a'e*'^ sinh^ (58) 

S\g)o/ S{g) — (T cosh^ — ae~^^ sinh^ (59) 

D\a)aD(a) =a + a (60) 

D\a)a'^ D{ol) =a'^+a*. (61) 



(49) We obtain 



M?\d4>^'^:rl)l = ^.,m<\,,,A.,,0 - km-lV^X 



X B 



(0* 



a,^,^^,Az,r7,(/> + ^i,m+lVm + l^a,<e,^z,Az ,r7,(/ 



(62) 

where the A J^^^^ ^^^ ^^ and ^J^^,^^,;,^,^,^ quantities are de- 
fined as 



^L:^,^z,Az,r,,0 = ^[(507O7r -tk^^tD] 



(63) 



Ai) 



.eff/ 



<f,,„A,.,> = cosher (5^70 - e'" sinhgfid^^n 



(64) 



Note that the €ij^rn,n coefficients present the following sym- 
metries. 



■m,n,m,n ^^ U 


(65) 


^i,j,m,n ^ ^m,j,i,n 


(66) 


^i,j,m,n ^ ^i,n,m,j 


(67) 



By inserting Eqs.(l50l)-(l62l) in Eq.(l48l) and by exploiting the 
symmetries of the eij^rn,n coefficients we obtain 

oo 
m,n=0 

(68) 
The QFI i^ampi(<^, 6>, (/), ^, ^, A, 7^) of the scheme is obtained 
by replacing gn -^ —g and gv -^ —g- This choice of the pa- 
rameters is equivalent to the case described in the main paper 
(with gH ^ —g, 9v ^ 9 and the additional 7r/2 phase shift 
in the probe state) leading to the same expression for the QFI. 
We finally obtain 



H{a,e,(t),£^,g,\r]) 



2\a\\r^ 



1 + 4/^(1 — rf) sinh g 



{ cosh[2{g - ^eff)] - cos(A + 2(^ - 26) sinh[2(^ - ^eff)]} 

(69) 
The optimal condition corresponds to the case cos (A + 2(/) — 
26>) = -1, where the QFI is 

g2(^-^eff) 

H^mA^.i.g.T]) = 2\a\^^r] ^ (70) 

1 + 4/^(1 — 77) sinh g 



Again, the dependence of the QFI H of (1691) on the parameter 
(j) to be estimated implies that to achieve its maximum i^ampi, 
an adaptive strategy (see Sec.|IIl| is necessary. 



VI. CLASSICAL FISHER INFORMATION FOR THE 
PHOTON-COUNTING MEASUREMENT 

In this section we describe the calculation for the clas- 
sical Fisher information associated with our scheme when 
photon-counting measurements are performed [Fig 121. The 
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Amplification Lossy photon- 
output state counting detection 

FIG. 2. Scheme of the different steps of the amplifier-based phase 
estimation protocol when photon-counting measurements are per- 
formed. 

output state of the protocol is described by the density matrix 
/^0 ^' while the measurement operators that describe photon- 
counting detectors are the projectors over Fock states 

nnw,n(-)=<fi)^<Tv) 



where li^^^,^ = |n^^))n(n^^^ |, with I = H,V labeling the op- 
tical mode. The probability distribution of the measurement 
outcomes can be evaluated as 

p(n(^),n(^)|0) = TV[p^'«'''n„,«,,„,.)] (72) 

The classical Fisher information associated to the probabil- 
ity distributions of the measurement outcomes is given by the 
following expression ||2|] : 



i^= E 



n,m=0 



p(nm,n(^)|(/)) 



[90P(' 



(73) 



For the amplifier-based protocol, the probability distribution 
p{n^-^^ , n^^^\(})) can be separated in two independent single- 
mode contributions as 

p{n^^\n^^^\(fy = \[ p(n(^)|(/)) (74) 

Here, p^^^ are the single-mode density matrices for modes / = 
H, V and 



p{n^'^\ci^) = Tr[p^'^Ii%] 



(75) 



In this case, the classical Fisher information can be separated 
in two single-mode contributions 



l=H,V 



(0 



(76) 



where 



«U^12 






(77) 



(71) 



A. Photon-number distribution of the amplified coherent 
states 

We begin by calculating the photon-number distribution of 
the amplified coherent states. The density matrix of the output 
state before the measurement stage is given by 

P^'^'^ = DH{lH)Dv{jv)SH{gf)Sv{gf) [p'^{Nf)^ 

(78) 
to evaluate the photon-number distribution, we exploit the fol- 
lowing identity between the elements of the density matrix ex- 
pressed in the Fock basis p = X]^^=o Pn,m\n){m\ and the 
Wigner function of a general single-mode state p, 

/CO />oo 
/ dxdpWp{x,p)Wn,m{x,p) (79) 
-cx) J —00 

where Wn^m (^, p) is the Wigner function associated to the op- 
erator \n){m\. Here, the (x, p) operators are defined according 



to A^xA^p > 1/16. The corresponding photon-number dis- 
tribution can be recovered from the diagonal elements pn,n^ 
by exploiting the expression of the Wigner function of a Fock 
state: 

WnA^^P) = -(-l)"i„[4(x2 +p2)]e-2(.^+P=) (80) 

TT 

Since the density matrix of the state p^ '^ = p\l ^ p\k is 
separable between the two modes, we can evaluate the distri- 
butions for the two components p\^ separately. The first step 
is the evaluation of the Wigner function for the single-mode 
density matrix: 



SI) _ 



effx ^/i/7vreff\Qt/ effN nt 



Di{^i)Si{gr)pT{Nr)Sj{gf)D]{^i) (81) 

The Wigner function for this state takes the form of a Gaussian 
distribution of the form 

,rr , X 2 1 -J-ij^['2(xi-^1)(pi-p1)<y7\ 



X e 



-^^^[(:.z-x?)Vr+bz-rf)vn 



(82) 
where the first order and the second order moments are, re- 
spectively 



X? = Re[7d 
P? = Im[7^] 



(83) 
(84) 



and 



(jf ^ = cosh(2^f^) + cos A^ sinh(2^f^) (85) 

(jf = cosh{2gf) - cos Xi sinh(2^f^) (86) 

afP = smXismh{2gf) (87) 

Here, gf^ and Xi are respectively the absolute values and the 
phase of the squeezing parameters gf^. We can now proceed 
with the calculation of the single-mode photon-number distri- 
bution p{n^^'^ 1^), which can be evaluated from the integral 



p{n 



(Oi 



OO POO 



CX) ^ — OO 



dxidpiWp(i){xi,pi)Wn,m{xi,Pl) 

(88) 



We first begin by performing the following rotation on the 
quadrature variables {xi^pi) -^ i^'uPi) of the Wp(i){x^p) 
function: 

x[ = xi cos iIji + Pi sin i/ji (89) 

p[ = -xi sin tlji + pi cos i/ji (90) 

x[^ =x[^ cos V^z + pj ° sin V^^ (91) 

p[ ° = -xj sin V^z + pj ° cos ^i (92) 

where ?/^/ = Xi/2. The Wigner function in this rotated quadra- 
ture set is 



w // /^ 2 1 -^^[(.;-.;0)^e^^r ] 

W^(i)(xi,p)) = — -™e ^+^"^1 



(93) 



X e ^ 



^^[(p;-pr)^e-^r"] 



The same rotation is performed on the Wn^ni^hPi), which 
presents radial symmetry and hence its form is not affected by 
the rotation according to 



TT 

(94) 
We can then proceed with the evaluation of the integral ([88l) . 
By performing the basis rotation {xi^pi) -^ i^'uPi) in the 
integration variable we obtain 



p{n^^^\^) =7r 



CO />oo 



dx[dpiW ii){x[,Pi)Wn,m{x[,Pi) 



CX) J — CX) 



(95) 

By expanding the Laguerre polynomials of the Wn,n{x[jp[) 
function we obtain 



p{n^^^\(j)) 



4(_l)n ^ 



E 



E 



(-4)'= fk 



/CX) />CX) 

-CX) J — OO 

X 6 ^ 

(96) 
The integrals in dx[ and dpj can be evaluated separately. We 
now define the following auxiliary functions 



A., = 1 



-B.r 



-^gf 



1 + 2Nf 



Cx, 



1 + 2Nf + e-29r 

{x'l^fe-^sf 
1 + 2iVfff + e-29r" 



^P' = 1 + 1 + 2Nf 
-B77, = 



1 + 2Arfff + e^sr 



^'"l + 2iVfff + e29r^^ 



(97) 

(98) 

(99) 

(100) 

(101) 

(102) 



where the B and the C terms depend on the phase (/>. Finally, 
by exploiting the definition of the confluent hypergeometric 
functions U{a, b; z), the single-mode photon number distribu- 
tion can be written as: 



2 in\ Ik 



Pinm = T^^e-^^-"^-^EE fc, ykl \j 



^ ' "^'^ k=Oj=0 

U[-j, 1/2, -2AAB,,f]U[-k + j, 1/2, -2ip,(4,)2] 



(A.,)^-+V2(A^jfc-.+i/2 



(103) 



B. Derivative of the photon-number distribution and classical 
Fisher information 

In order to evaluate the classical Fisher information accord- 
ing to Eqs. (17 611771) , we now need to evaluate the derivative of 
the photon-number distribution p(n'^^^ |0). The latter can be 
written in the following form 

n k 

^^ U[-j, 1/2, -2A,,{B,,f]U[-k^j, 1/2, -2^(^pJ'] 



VII. MODELING THE EXPERIMENT 

Here we discuss the theoretical model for the analysis of the 
experimental data of the protocol. In the implementation de- 
scribed in the main paper, no phase stabilization is performed 
on the optical path of the pump beam, hence the phase varies 
randomly at each experimental run. To model such effect, an 
average on the phase A with a uniform distribution V{X) = ^ 
must be performed on both the signal and the fluctuations. In 
this case, the average signal in the two polarizations H and V 
is given by 



{A,,y^y^{A,,)^-^^y' 



(104) 



Here, ujn,kj includes all the coefficients independent from the 
phase (j). The derivative of the photon-number distribution 
p{n^^'^ \(j)) can then be written as the sum of three terms 



{uh) = V [n + |ap^(l + 2n)cos2((^/2)] (112) 
J^ = r][n^ |ap^(l + 2n) sin2(0/2)] (113) 
The average number of the count rates (D) is then given by 



(D) = |a|^7/^cos^(l + 2n) 



(114) 



dMn^'^\^) = Y,DMn 



(Oi 



(105) 



The term Dpnin^^^ |^) presents the derivative of the exponen- 
tial e~'^^^''i^^pi\ leading to: 

DpH{n^'^^) = (-2)S^(C,, +CpJp(n^^)|(/)) (106) 

The terms Dpy{n^^^\(j)) and I)p3(n'^^^|^) exploit the follow- 
ing relation involving the derivatives of the confluent hyper- 
geometric functions: 

d^U[a, 6, /((/))] = -aU[a + 1,6 + 1, /((/>)] 5^ /((/>) (107) 

The remaining two terms can the be written as: 

I k=0 j=0 ^ ^ 



In the high losses regime investigated throughout the paper, 
the number of the photons effectively impinging on the detec- 
tor is smaller than one, since rj{n±) < 1. In this regime, the 
single-photon counting process is described by a Poissonian 
statistics. Hence, the fluctuation on the difference signal can 
be evaluated as 



cj\{D)) = a'iiriH)) + cr'iinv)) = W) + (ny) (115) 



By explicitly substituting the expressions for (uh) and (ny) 
we obtain the following expression for the sensitivity 

^^ _ |apev^(l + 2n)|sin0| 
'"P V2n+|a|2e(l + 2n) 



(116) 



The optimal point is achieved for ( 
tivity is 

|ap^V^(l + 2n) 



7r/2, where the sensi- 



^, J7[l-j, 3/2, -21., (E,,)2]j7[_fc + j; 1/2,-24^ (4J2] 
xi(-4)i,,B,,(a^B,,) 



S, 



exp 



V2n+|a|2^(l + 2n) 



(117) 



(108) 



and: 



n \ 2'= fk 



^^ U[-j, 1/2, -2A,,{B,,)']U[1 -k^j, 3/2, -2Ap,{Bp,f 

x{k-j){-4)Ap,Bp,{d^Bp,) 

(109) 
Finally, the classical Fisher information can be evaluated ac- 
cording to: 



Appendix A: Quantum Fisher Information 

Here we briefly review the properties of the quantum Fisher 
information for mixed states. Let us consider a family of states 
0-0 depending on a parameter (/). Such family of states can be 
exploited to estimate the value of the parameter (j). In local es- 
timation theory, the maximum amount of information that can 
be extracted on the parameter (j) with M repeated measure- 
ments is given by the QFI H^p. More specifically, the vari- 
ance of any estimator of the parameter (j) satisfies the quantum 
Cramer-Rao inequality: 



l=H,V 



S^^> 



1 



where: 



^r-E 



p(n(')|</)) 



(110) 



(111) 



MH^ 



(Al) 



Here, iJ^ represents the optimization of the classical Fisher 
information over all possible choice of the quantum measure- 
ment. In general, the quantum Fisher information of the fam- 
ily of states (jfj, is given by the following definition: 



n=0 



H^=Tr[a^Ll] 



(A2) 



where L^ is the symmetric logarithmic derivative of a^: 



Appendix B: Mathematical relations 



L(f)Cr(j) + (701/0 



(A3) 



By expressing the density matrix in terms of its spectral de- 
composition (70 = ^^ (7m|Cr7i)(Cm|, the quantum Fisher in- 

si): 



formation can be evaluated as I 



^^ = E 



2^en,m|(Cm|S0Cn)|' (A4) 



Here, 9^(7^ is the derivative of the eigenvalues with respect to 
(j), and \d(f)C,n) is the derivative of the eigenvectors written in a 
^-independent basis {\k)}'. 



I^^Cm) = ^{d^Cmk)\k) 



(A5) 



Finally, the coefficient en,m is given by the following expres- 
sion: 



((7n - (7mY 



(A6) 



In this Appendix we report some mathematical relations ex- 
ploited in the calculation of the Fisher information. 

Thermal state. - The thermal single-mode state is defined 
as: 

-J CX) 

p'n7V) = — =^x"|n)(n| (Bl) 

n—0 

with X = A/'/(l + N), where N is the average number of 
photons of the state. 

Lossy squeezed vacuum. - The state generated by the ac- 
tion of a lossy channel on the squeezed vacuum state can be 
written as according to Ij] : 

Cr, [S{gm{0\SHg)] = S\g''')p'\N''')Si9''') (B2) 

The effective modulus of the squeezing parameter g^^ and the 
effective thermal noise N^^ take the form: 



^'"^ = 7^47) 



P 



N' 



■eff 



-1 + /PM 



where: 



P = r]e^^ + 1-7/ 
M = r]e~'^^ + 1-7? 



(B3) 
(B4) 



(B5) 
(B6) 
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